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Abstract 

The Helmholtz equation for symmetric, traceless, second-rank tensor fields in three-dimen- 
sional flat space is solved in spherical and cylindrical coordinates by separation of variables 
making use of the corresponding spin-weighted harmonics. It is shown that any symmetric, 
traceless, divergenceless second-rank tensor field that satisfies the Helmholtz equation can be 
expressed in terms of two scalar potentials that satisfy the Helmholtz equation. Two such 
expressions are given, which are adapted to the spherical or cylindrical coordinates. The 
application to the linearized Einstein theory is discussed. 
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1 Introduction 



The gravitational field in the linearized Einstein theory can be represented by two symmetric, 
traceless, second-rank tensor fields which, assuming that the fields vary harmonically in time, 
obey the Helmholtz equation outside the sources. As in the case of other nonscalar equations, 
the solution of the Helmholtz equation for second-rank tensor fields in noncartesian coordinates 
is a difficult problem owing to the coupling of the fields components. However, when a nonscalar 
equation is written in spherical or cylindrical coordinates, a considerable simplification can be 
obtained by using spin- weighted quantities and the spin- weighted harmonics. 

In this paper the Helmholtz equation for spin-2 fields (i.e., symmetric, traceless, second- 
rank tensor fields) is solved by separation of variables in spherical and cylindrical coordinates 
making use of the spin- weighted harmonics. A similar treatment for the vector (spin-1) Helmholtz 
equation is given in Refs. PJEHH]- In Sect. 2 the Helmholtz equation for spin-2 fields is solved 
in spherical coordinates and it is shown that the divergenceless solutions of this equation can be 
expressed in terms of two scalar (Debye) potentials that satisfy the Helmholtz equation. The 
expressions for the divergenceless solutions in terms of potentials obtained here are equivalent to 
those found in Refs. [UO2] for the multipoles with j > 1. In Sect. 3 the Helmholtz equation for 
spin-2 fields is solved in cylindrical coordinates and an expression for the divergenceless solutions 
in terms of potentials adapted to the cylindrical coordinates is obtained. In Sect. 4 two alternative 
expressions for the solutions of the linearized Einstein vacuum field equations in terms of scalar 
potentials are given. 



2 Separation of variables in spherical coordinates 

Let {ei, e2, 63} be an orthonormal basis, a quantity ij has spin- weight s if under the transformation 

ei + ie 2 -> e ta (ei + ie 2 ) (1) 

it transforms according to 

77 -► e isa r). (2) 

The five independent components of a symmetric, traceless tensor, tij, can be combined to form 
the quantities 

t±2 = i(tll-t22±2ttl 2 ) = i(t33 + 2*ll±2itl2) ) 

t±i = =Fj fas ± itas) , (3) 
1 

*o = 2* 33 ' 

so that t s has spin-weight s. If the components tij are real, then 

t s = (-l) s t_ s , (4) 

where the bar denotes complex conjugation. Similarly, in the case of a vector field F, the 
combinations 

F±i = (F ± %F 2 ) , F = "^F 3 , (5) 
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have spin-weight ±1 and 0. 

By choosing the basis {ei, e 2 , 63} as the basis {eg, e^, e r } induced by the spherical coordinates, 
one finds that the Helmholtz equation for a symmetric, traceless tensor field t, 



V 2 * + k 2 t = 0, 



(6) 



written in terms of the spin-weighted components © amounts to 



^d r (r 2 d r t +2 ) + ^St + 2 + ^t +1 + k 2 t + 2 = 0, 

\d r (r 2 d r t + i) - -^t+i + \ddt + i - 4rSi+2 + 4rS*o + k 2 t + i = 0, 
1 6 12 

-^d r (T 2 d r t^)- -^t + ^3fit + -5 (St-i - + A: 2 t = 0, 

-^a r (r 2 a r t_i)-^t_i + -^gg*_i + ^gt_ 2 -^§*o + ^ 2 *-i = o, 

\d r (r 2 d r t- 2 ) + - + ^-2 = 0, 



(7) 



where, acting on a quantity 77 with spin- weight s |3[7], 



97? = -sin s 6> [dg + -— 9,4 Wry sin s 9), 
s sin ft 



0?/ 



sin" 5 e a 



sml 



(r)sin s 9). 



(8) 



(The expressions Q can be readily obtained by using the spinor formalism of Ref. [H].) 
We seek separable solutions of Eqs. (JZJ) of the form 



t±2 



j(j + 1) 



1/2 



.(i-i)0' + 2) 

*±i = [?(3 + l)] 1/2 ^±i(r)±i^m(e,0), 
*o = j(j + l)go(r)Y jm {6,(f)), 



(9) 



where j is an integer greater than 1 and s Yj m are spin- weighted spherical harmonics [|3E|- The 
tensor field given by Eqs. © is an eigentensor of J 2 and J3 with eigenvalues j(j + 1) and m, 
respectively (see, e.g., Ref. P). Substituting Eqs. © into Eqs. (J7| one obtains the set of 
ordinary differential equations 



d 2 | 2 d (j-l)(j + 2) | fc2 

dr 2 r dr r 2 

^ 2 | 2 d j(j + l)+4 | fe2 

dr 2 r dr r 2 



4(j-l)(j + 2) 

g±2 h ^3 s^ 1 

1 3j(j + 1) 

9±i + -^5±2 H 5 5o 



d 2 , 2 d j(j + 1) + 6 fc2 



dr 2 ^ r dr 



2 , 

5o + -2(0-1 +9+i, 



0, 
0, 
0. 



(10) 
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By combining Eqs. (jTU|) we find decoupled equations for [H] 

g +2 -g-2 - 2(j + 2)( 5+ i -9-1), 

g +2 -g-2 + 2(j -l)(g +1 -g-l), 
g+2+g-2 - 4(j + 2)( 9+1 + + 6(i + l)(j + 2) 50 , 
g +2 +g-2 - 2(g +1 +g-l)-2j(j + l)g , 
5+2 + 5 -2 + 4(j - 1)( 5+1 + 9-1) + 6j(i- 1)50, 
whose solutions are spherical Bessel functions provided k ^ 0. Thus, from Eqs. © we get 

t±2 = \ [{j ~ + 1)0' + 2 )] 1/2 { aj j+ 2(kr) + bn j+2 (kr) - 2[cj j (kr) + dn^fcr)] 
+ejj- 2 (kr) + fnj- 2 (kr) ± 2 [-Aj j+1 (kr) - Bn j+1 (kr) 
+Cjj-i(kr) + Dnj-i(fe-)]} ±2^m, 

t±i = - + 1)] 1/2 { -(i + 2) [ai i+2 (A:r) + bn j+2 (kr)} + c jj(fcr) + dn^kr) 

+(j - 1) [ej^ 2 (kr) + /n J -_ 2 (fcr)] ± (j + 2) L4^+i(Ax) +Bn J+ i(ir)] (11) 
±(j - 1) [Cjj-^kr) + Dn j ^ 1 (kr)]} ±1 Y jm , 

to = | + 1} 2 (J + 2) [a jj+2(fcr) + bn j+2 (kr)} + [ CJj (fcr) + dn^kr)] 

+ [eii- 2 (A:r) + f nj -2(kr)]\ Y jm , 

where a, 6, c, d, e, /, ^4, B, C, and -D are arbitrary constants. 

The cases j = 1 and j = must be treated separately since s Yj m = for |s| > j. We find 
that, also in these cases, the separable solutions of Eqs. (J7J) are given by Eqs. Qll|). 

As in the case of the vector Helmholtz equation, the fact that the radial equations can be 
decoupled is related with the existence of an operator that commutes with J 2 , J3 and V 2 + k 2 
[2]. Such an operator can be chosen as 

K = L k S k + 2, (12) 

where L k and S k are the operators corresponding to the cartesian components of the orbital and 
spin angular momentum, respectively, and the summation convention applies. 

For a symmetric, traceless tensor field t, the spin- weighted components of Kt are given by 

(Kt) +2 = 2g* +1 , 

(Kt)+i = --Si+ 2 -3t+i + ~Si , 

(Kt) = -6i+i - 4t + 3t_i, (13) 

(Kt)-i = ^St_2 - 3t_i - ^m , 

(Kt)- 2 = -29i_i. 

Using Eqs. Q (which give the eigentensors of J 2 and J3) and (j!3[) one can readily find the 
common eigentensors of J 2 , J3 and K; for these tensor fields only one of the five radial functions 
g s (r) is independent and, therefore, Eqs. (|T0|) reduce to a single equation (cf. Ref. |2j). 
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Since J 2 = L 2 + 2LkSk + S 2 and, for spin-2 fields, S 2 t = 6i, it follows that 

K = \(J 2 -L 2 -2\ (14) 

which shows that the eigentensors of J 2 ,5 2 , and K are also eigentensors of L 2 and using Eqs. 
(jl3H14|) it is easy to see that the separable solution (jllj) is a superposition of five eigentensors of 
L 2 with eigenvalues 1(1 + 1), where I coincides with the index of the spherical Bessel functions 
appearing in Eqs. Ijlljl . Assuming that under the parity transformation, r — > — r, e r and 
are left unchanged and eg changes sign and taking into account that s Yj m is transformed into 
(— 1) J - s Yj m , one finds that the five eigentensors of L 2 contained in Eqs. (jllj) are also eigentensors 
of the parity operator with eigenvalue (—1)'. 

The divergence of a second-rank, symmetric, traceless tensor field, t, is the vector field div t 
whose cartesian components are given by (divt)j = djtij, where dj = d/dxj. The components of 
the divergence of t with respect to the basis {eg, e^, e r } are determined by [see, e.g., Ref. jSj, Eq. 
(44)] 

(div t) s = i= {--&.+! ~ ^d r (rh s ) + -gt s _i) , (15) 

with the spin- weighted components of div t defined as in Eq. (J5J). Substituting Eqs. (fTTj) into 
Eqs. Ijl5j) and using the recurrence relations for the spin- weighted spherical harmonics and for 
the spherical Bessel functions, one finds that the separable solution of the Helmholtz equation 
given by Eqs. Ijlljl has vanishing divergence if and only if 

j(2j - l)c , _ j(2j -l)d „ _ (j + l)(2j + 3)c 



3(j + 2)(2j + l)' 3(j + 2)(2j + l)' 3(j - l)(2j + 1) ' 

= 0; + l)(2 J + 3)d A = J^1 B = i^l D . 

3(j-l)(2j + l)' j + 2 j + 2 

Substituting Eqs. Ijl6j) into Eqs. Ijlljl and making use of the recurrence relations for the Bessel 
functions one gets 



-^cV^i + -( \d 2 r 2 - k 2 ) CISV2 
r 2 \ r 



ik - 1 
t+\ = — SSS^i - — -^d r rddd^2, 
2r 2r z 

t = -Lm&h, (17) 
2r z 

ik 1 
t-i = — SSS^i + — Td r rddd^ 2 , 
2r 2r z 



where 



^ = % M^l C ^ kr )+ Dn j( kr )] Y jrn, (18) 



^ 2 55 B(j l c Mkr) + d nj (kr)] Y jm . 
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Clearly, the functions ipi and ip2 are separable solutions of the scalar Helmholtz equation. It may 
be noticed that Eqs. (|17|) contain no reference to the value of j. On the other hand, for j = 1, 0, 
from Eqs. and (jl5[) one finds that if k ^ and div t = then, necessarily, t = 0. 

Thus, by virtue of the completeness of the spin- weighted spherical harmonics, any divergence- 
less solution of the spin-2 Helmholtz equation © can be expressed as a superposition of separable 
solutions of the form (|17|) where, now, ipi and ip2 are two solutions of the scalar Helmholtz equa- 
tion that are superpositions of solutions of the form (|18|) . In view of Eq. if ipi and ip2 are 
real then t is real. (The factor i was included in Eqs. (|17|) in order to produce this relation.) 

The components (|17jl can be written in terms of certain tensor operators Uij,Vij [3], whose 
cartesian components are defined by 

Uij = LiXjip + LjXiip, Vij = e imn d m U nj (ip), (19) 

where now 

L = rxV, X = VxL-V. (20) 

It is easy to see that for any well-behaved function ifi, Uij(tjj) and Vij(^) are symmetric, traceless, 
divergenceless tensor fields. By computing the spin- weighted components of Uij(ip) and Vij(ip) 
with respect to the basis {ee,e^,e r } [S] one finds that the expressions ()17j) are equivalent to 

Uj = kU ij (ih) + Vi j (fh)- (21) 

Equations (fT9*|) imply that 

BimndmVnM) = -^i(V^), (22) 

and therefore the tensor field ()21|) satisfies 

(curl t)ij = e imn d m t nj = k 2 Uij(ip 2 ) + kVij(i>i). (23) 

According to Eqs. (|17ft . the scalar potentials generating a divergenceless solution of the spin-2 
Helmholtz equation are determined by 

8695^2 = 2r 2 i , (24) 

and, by comparing Eqs. (|21j) and (|2*3*)) . 

2r 2 

gggg^i = ^(curi tw (25) 

The uselfulness of Eqs. (|24|) and (|25[1 can be illustrated by obtaining the expansion of a 
circularly polarized plane wave in spherical waves. The cartesian components of a spin-2 field 
corresponding to a circularly polarized plane wave with helicity ± propagating in the z-direction 
are given by 

/ 1 ±i 

it;/) A ±i -1 ) e ikz , (20! 

V 00 
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where A is a constant. From Eqs. © and (|2fi|) one finds that 2r 2 to = XiXjtij = A(r sin 9 e^ 1 ^) 2 e tkz ; 
therefore, by using the expansion of e tkz in terms of spherical harmonics and the recurrence 
relations for the associated Legendre functions (see, e.g., Ref. ^01), Eq. ([21]) yields 

9999^2 = E [ 47r ^' + ^O' " + !)0' + 2 )] 1/2 Ji(^) 
hence, we can choose 

1/2 



: ,±2, 



(27) 



02 



1 



i=2 



4vr(2j + 1) 



v> jj(kr)Y j}±2 . 



(28) 



.(i-i)i(i + i)(i + 2). 

Since the tensor field ()26j) is such that that curl t = ±kt, from Eqs. (|24[1 and (|25[1 we see that 

V>i = ±^2. (29) 
By substituting Eqs. (|2"8)) and (|2"9*|) into Eqs. ((T7|) or (|2"Tj) one gets the desired expansion (cf. Ref. 
™- 

In the case where k = 0, the Helmholtz equation reduces to the Laplace equation and by 
assuming a separable solution of the form Q we obtain [cf. Eqs. pip] 

\ [(J " l)i(i + + 2)] 1/2 { o^ 2 + for--''- 3 - 2[cr J ' + dr^" 1 ] 
+er j ~ 2 + fr- j+1 ± 2 [-Ar j+1 -B r - j - 2 + Cr j ^ +Dr~ j } } ± 2 Y jm , 

~ + 1)] V2 {-U + 2 ) [° rj+2 + &^ j_3 ] +cr j + dr-^ 1 
+ (j - 1) [er j - 2 + fr- j+1 ] ± (j + 2) [.4r J ' +1 +Br^ 2 ] 
±(j-l) [Cr^ + Dr^]} ±1 Y jm , 



t±2 



t±l 



(30) 



to 



(J + 1)(J + 2) r +2 3] j(j + 1) r , 
I a r J ^ + o r J J H [cr + a j 



+ 



i(j - 1) 



e r 



j-2 



+ fr 



-i-il 



for j = 0,1,2,... 

Substituting expressions (j30J) into Eqs. p5j) one finds that the solution of the Laplace equation 
given by Eqs. (j30j) has vanishing divergence if and only if 

a = c = d = f = A = D = 0. (31) 

When these relations are inserted into Eqs. (|30j) . for j > 1, they can be written as [cf. Eqs. p7|l ] 

t+2 = 



1 



t+i 

to 

t-1 

t-2 



398-01 
9999-02, 



1 

272 



<9 r r 99902, 



2r 
1 

2T 2 " 
2r C 

-=fl r r 2 93i/>i + 7; — 7:d 2 v 2 890 2 , 
r z 2r 



(32) 



SSV'i + 7-^5^99902, 
2r z 
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where 



tp2 



(7 -1)0' + 2) 
1 

0-i)0 + 2) 



[(j-l)Cr^(j+2)Br-H] 



(33) 



J-l 
3 + 1 



er + 



J + 2 



Y, 



In this case, ipi and ■02 are separable solutions of the scalar Laplace equation. 

From Eqs. (|3UH31|) we see that the divergenceless solutions of the Laplace equation with j = 1 
are given by 

r t±2 = o, 

3 = 1 : t ± i = ^ (±Br 
I to = 3br~ 4 Y lm , 
which can be written in the form (|32|) with 

iB . 

- — (>Wnr + Sh r , 



>r~ 3 -br 4 ) ±iY lm , 



(34) 



where h r 



is a solution of 



ipl = 

4>2 = 



b 



r (Yi m lnr + 3h 

n 



Y lr 



(m = ±1,0). 



(35) 



(36) 



Owing to Eq. (|3l))) . the scalar potentials (|33|) satisfy the Laplace equation. Finally, in the case 
where j = 0. Eqs. (130131 1) yield 



j = 

This solution is of the form (|32|) with 

^i = 0, 




V'2 



0, 

o, 

br- 3 Y 00 . 



Yqq In (r esc I 



(37) 



(38) 



which are solutions of the Laplace equation. It may be noticed that the scalar potentials (|35|) 
and (|38jl diverge at 6 = 0, tt, and are not separable. 

Thus, any divergenceless solution of the spin-2 Laplace equation can be expressed in the form 
(|32|) . where and ip2 are solutions of the scalar Laplace equation. Equations (|32p are equivalent 
to 



(39) 



3 Separation of variables in cylindrical coordinates 

Taking now the basis {ei,e2,e3} as the orthonormal basis {e p , e0,e 2 } induced by the circular 
cylindrical coordinates, the Helmholtz equation for a symmetric, traceless, second-rank tensor 
field is equivalent to the set of uncoupled equations 

d 2 z t s + mt s + k 2 t s = 0, (s = ±2,±l,0), (40) 
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where acting on a quantity rj with spin- weight s [2], 

Sr? = ,r (d p ■ -0, ; 



P 

Brj = -p-fdp-fyofr,). (41) 

We seek solutions of Eqs. 1)40(1 of the form 

t s = 9s(z) s F am (p,^), (a = ±2, ±1,0), (42) 

where the s F am are spin- weighted cylindrical harmonics 3 . A tensor field of the form 1)42(1 is 
eigentensor of J3 and of the square of the linear momentum perpendicular to the z-axis, P? + P 2 , 
with eigenvalues m and a 2 , respectively Substituting Eqs. ((4^|) into Eqs. (jlU]) one finds 

^-7 2 )^ = 0, (* = ±2,±1,0), (43) 
where 

7 2 = a 2 - k 2 , (44) 

therefore, if 7 7^ 0, g s (z) = A s e lz + B s e~ JZ and if 7 = 0, g s (z) = A s + B s z, where A s and B s 
are arbitrary constants. Thus, assuming that a is different from zero, Eqs. 1)40(1 admit separable 
solutions of the form 

U = (A s e^ z + B s e^ z )[C s ( s J am ) + D s ( s N arn )}, (7 + 0), (45) 

and 

t s = {A s + B s z)[C s ( s J am ) + D s ( s N am )], (7 = 0), (46) 
where A s , B s , C s and .D,, are arbitrary constants and [3] 

s Jam(p, 0) = J m+S {ap)e im * , s iV Qm (p, 0) = N m+S {ap)e im * , (47) 

with and N u being Bessel functions. For a = 0, the functions s F am diverge when p goes to 
zero or to infinity, or they do not vanish when p goes to infinity. 

The components of the divergence of a symmetric, traceless, second-rank tensor field with 
respect to the basis {e p ,e^,e 2 } are given by [see, e.g., Ref. |Hj, Eq. (44)] 

(divt), = {-Bt s+1 - 28 z t s + , (48) 

therefore, the tensor field 1)45(1 has vanishing divergence if and only if 



§ (A C + A ±2 C ±2 ) = jA ±1 C ±1 , %(A D + A ±2 D ±2 ) = 7 A ±1 D ±1 , 

q(B C + B ±2 C± 2 ) = -jB ±1 C ±1 , %(B D + B ±2 D ±2 ) = -jB ±1 D ±ll 

f (A +1 C +1 + A_xC_i) = 7 A C , %{A +l D +1 + A_iZ>_i) = 7A0A), 

f (B+iC+i + = - 7 £oCo, f (S+iD+i + = -"/B D . 



(49) 
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Introducing the combinations 

01 = ^(A+iC+i - A_iC_i), a 2 = ±{A +1 D +l - A^D^), 

h = ±(B +1 C +1 - B_!CLi), 6 2 = ^(B+i^+i - B-1B-1), 

a 3 = 2^(A +1 C +1 +^_ 1 C_ 1 ), a 4 = ^(A+xD+a + ^LU), ^ 

^3 = -3^(B+iC + i + B-1C-1), 6 4 = -^(B+iD+i + B-iD-i), 



and assuming that the conditions (|49[) hold, one finds that the components (|45f) can be written 
as [cf. Eqs. (Q2J) and 

i - 1 - 

i+i = -999^1 - -d z dddif) 2 , 

t = ^9999^2, (51) 

% 1 

f_i = -999-01 + 2^999^2, 

t_ 2 = ^99^1 + ^(^-^)99^2, 



where 



1P1 = ^[(aiey z + b 1 e-^) J am + (a 2 e'y z + b2e-^) N arn ], (52) 
2 

V>2 = ^ [( a 3 e ^ + &3e~ 7Z ) O^m + (a4e 7 ^ + M" 7 *) O^am] , 

which are solutions of the scalar Helmholtz equation. 

Similarly, one finds that if the field given by Eqs. (|46|) has vanishing divergence then its 
components can be written in the form (|51|) . where 

2i 

Ipl = —n [(Ol + hz) J am + (o 2 + &2^) O^om] , (53) 

2 

"02 = -n [(a 3 + 6 3 ^) Jam + (a4 + M) O^om] , 

cr 

which satisfy the scalar Helmholtz equation, and 



ai = ^(A+iC+i - A_iC_i), a 2 = ^(A+iD+i - 
6i = ^(B +1 C +1 - B_ 1 C_ 1 ), 6 2 = 2k(B +lJ D +1 - B-iD-x), 
as = -0^(^+2(^+2 + A_ 2 C_ 2 ), a 4 = -^(A +2 D +2 + A-2B-2), 



(54) 



63 = ^(A +1 C+i + A_iC_i), 64 = ^(A +1 D +1 + 
It can be shown that Eqs. ()51() are equivalent to 

tij = Wi j (rp 1 )+Z ij (7p 2 ), (55) 
where the tensor operators Wij and Zjj are given in cartesian coordinates by 

Wij (V) = MiNj iJj + MjNi 1P, Zij (V) = e imn d m W nj (V) , (56) 
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with 

M = e z x V, N = V x M, (57) 

[cf. Eqs. (|19I21() ]. It is easy to see that for any well-behaved function ip, Wij(ip) and Zij(ip) are 
symmetric, traceless, divergenceless tensor fields, and that 

EimndmZnM) = -Wy(VV)- (58) 

The solutions (|42|) with a = can also be written in the form (|51|) . in terms of two scalar 
potentials ipi and ip2 that satisfy the Helmholtz equation but they are not separable. Owing to 
the completeness of the spin-weighted cylindrical harmonics, any divergenceless solution of the 
spin-2 Helmholtz equation can be written in the form (|51j) . where ipi and ip2 are solutions of the 
scalar Helmholtz equation. If ip\ and ip2 are real, then tij is real. 



4 Application to the linearized Einstein equations 

The Einstein vacuum field equations linearized about the Minkowski metric can be written in 
cartesian coordinates in the gauge-invariant form 



diEij = 0, d l B i3 ={), (59) 

19 Id 

~C~dt — ^-imnOm^nj i ~c~~3t* 



I — Ei mn d m B n j, pi + Bij — Ei mn d m E n j , (60) 



where Eij and B^ are symmetric traceless tensor fields defined by 

Eij = K()iQj , Bij = ~Koimn Ejrnm (61) 

and K a p~ f $ is the curvature tensor to first order in the metric perturbation (see e.g., Refs. |Hll2j). 
Prom Eqs. (fBTH IfiU]) it follows that the fields Eij and Bij obey the wave equation; therefore, 
assuming that Eij and B^ have a time dependence of the form e~ lujt , the fields Eij and B^ 
satisfy the Helmholtz equation with k = tu/c. According to the results of Sect. 2, if u; 7^ 0, Eij 
can be expressed in the form 

E l j = kUij(^ 1 ) + V l j^ 2 ), (62) 

[cf. Eq. (|21|) ] where tpi and ip2 are regular solutions of the scalar Helmholtz equation. From Eqs. 
(|Hn|) and (|2"3*|) we see that the field B^ corresponding to (|6*2"|) is given by 

B^ = -i[kUij(^ 2 ) + Viji^i)]. (63) 

In the static case (u = 0), E^ and Bij must satisfy the Laplace equation; hence, Eij can be 
expressed in the form Eij = Uij(ipi) J t-Vij(iJj2), where ipi and ^2 are solutions of the scalar Laplace 
equation [Eq. (|3H|) ]. Equations (|6T)|) and (|22|) give = £imnd m E n j = Vij (V'l), which implies that 
Uiji^x) = [see Eqs. (|32I3Sj) ]. Thus 

Eij = Vij{^2)- (64) 

In a similar manner, it follows that 

Bij = Vij{^), (65) 
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where ipn is a solution of the scalar Laplace equation. (It may be noticed that Eqs. Ijfi4lfi5jl can 
be obtained from Eqs. (|62W63j) by simply setting k = 0.) 

Alternatively, the components E^ and B^j can be expressed in the form (|55j) . Assuming again 
that the time dependence of the fields is given by a factor e _iW *, in the case where w^O, 

E ij = kW ij (il>i) + Z ij (ih), (66) 

where ifti and -02 are solutions of the scalar Helmholtz equation and the factor k has been intro- 
duced for convenience. Then, Eqs. (jfiU]) and (jSHj) imply that 

Bij = -i[kWij{$2) + Zijifo)]. (67) 

On the other hand, when to = 0, 

= Zij^t), B^ = Zij(ip 2 ), (68) 

where tpi and ip2 are solutions of the scalar Laplace equation. 

In the standard approach, the Einstein vacuum field equations linearized about the Minkowski 
metric are written in terms of the metric perturbations, which are affected by the gauge transfor- 
mations induced by the infinitesimal coordinate changes. By contrast, the curvature perturbations 
Kafirs j which are equivalent to the fields Eij and B^ , provide a gauge- invariant description of the 
gravitational field Eqs. (|62ll63j) . which are adapted to the spherical coordinates, yield a multipole 
expansion of the gravitational field. The gravitational potentials appearing in Eqs. ()62I68|) for 
fields generated by localized sources can be expressed in terms of the energy-momentum tensor 
of the sources as in Refs. 

5 Concluding remarks 

It is known that in flat space-time a massless field of arbitrary spin can be expressed in terms of two 
real potentials or of a single complex scalar potential (see, e.g., Ref. ^3]); the results presented 
above and in Refs. pj |H] show specifically that there exist operators such that when applied to 
a solution of the scalar Helmholtz equation which is eigenfunction of J 2 (= L 2 ) and Ja(= L3). 
they yield solutions of the spin-1 and spin-2 massless field equations that are eigenfunctions of J 2 
and J3 with the eigenvalues of the potential and that, similarly, there exist operators that map a 
solution of the scalar Helmholtz equation which is eigenfunction of P 2 + P| and J3 into solutions 
of the spin-1 and spin-2 massless field equations that are eigenfunctions of P 2 + P| and J3 with 
the same respective eigenvalues. 
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